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We discuss the use of a hysteretic Josephson junction to detect current fluc- 
tuations with frequencies below the plasma frequency of the junction. These 
adiabatic fluctuations are probed by switching measurements observing the 
noise- affected average rate of macroscopic quantum tunneling of the detector 
junction out of its zero-voltage state. In a proposed experimental scheme, 
frequencies of the noise are limited by an on-chip filtering circuit. The third 
cumulant of current fluctuations at the detector is related to an asymmetry 
of the switching rates. 
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1. INTRODUCTION 

Complete description of charge transport through a mesoscopic con- 
ductor requires knowledge of the probability distribution of current through 
the conductor. In general, such distributions are characterized not only by 
the average current and variance, but also by higher order moments. This 
study of fluctuations has attracted intense theoretical activity in the re- 
cent years, and powerful frameworks such as full counting statistics (FCS) 
have been formulated.'^ Starting from fundamental microscopic theories, 
moment-generating functions are known for several systems such as tunnel 
junctions and quantum dots. Experimental investigation of the fluctuations 
beyond the variance is, however, not so well-established. Experiments ob- 
serving higher moments of current or voltage have remained difficult and 
time consuming as filtering and bandwidth requirements are hard to ful- 
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fill — detection of higher moments of fluctuations is typically characterized 
by very weak signals and wide bandwidth measurements performed at low 
temperatures where thermal effects are suppressed. 

The first measurement^ of the third moment of current fluctuations 
across a voltage-biased tunnel junction support the theoretical concepts, al- 
though the use of conventional amplifiers and mixers requires long averaging 
times. Furthermore, the electrical environment of the conductor can signif- 
icantly affect the measured statistics.'^ In a more recent measurementj^ the 
use of a slightly different detection scheme provides a more direct access to 
fiuctuations of voltage up to the third moment. 

Among the growing number of experimental findings, remarkable results 
have been obtained using real-time detection of single electrons,'^ which can 
be utilized in the measurement of full counting statistics of electron trans- 
port. Recently, even further progress was accomplished in a measurement of 
the distribution of tunneling current through a quantum dot.'^ These inno- 
vative approaches are examples of on-chip detection of fiuctuations, which is 
the design goal of future experimental schemes as well. In this way, one can 
perform faster measurements on a wider bandwidth and avoid the problems 
with remotely connected amplifiers. 

On the other hand, the above techniques based on direct counting of 
electrons are best suited up to moderate frequencies and extremely low cur- 
rent levels in the pA range and below. Other experimental approaches are 
needed at considerably higher currents and for frequencies in the range of 
several GHz, indispensable for the characterization of many involved pro- 
cesses. One solution is the use of Josephson junctions (JJs) as noise probes. 
With small low-capacitance junctions, the sensitivity of the Coulomb block- 
ade can be used to characterize voltage fiuctuations generated by a meso- 
scopic conductor.'^'^ Another possibility is the use of larger current-biased 
junctions as threshold detector to probe the higher moments of current 
noise. In this case, switching rate of the junction from the supercurrent 
to the normal state depends strongly on the fiuctuations in the bias cur- 
rent. The applicability of such a JJ detector to measure shot noise has been 
demonstrated but convincing experimental results concerning the higher 
moments of fiuctuations have not yet been reported. 

A JJ appears to be quite an attractive threshold detector at the first 
view: it detects with a certain probability a current exceeding a given thresh- 
old, and switches consequently from the superconducting state to a well 
detectable normal state. A more detailed description of the dynamics of 
the junction acting as a noise detector reveals several mechanisms though: 
apart from the usual crossover from thermal activation (TA) to macroscopic 
quantum tunneling (MQT) as a function of temperature, the response of a 
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hysteretic Josephson junction to current noise depends inherently on both 
the spectrum and distribution of the fluctuations. In the TA regime, the 
switching of the junction occurs by thermally activated escape over a high 
potential barrier. On the contrary, at lower temperatures the escape results 
as the phase of the superconducting order parameter over the junction tun- 
nels quantum mechanically through the barrier. These processes are aff'ected 
by current fluctuations, and the influence of noise on the escape dynamics 
has been correspondingly analyzed for junctions in different parameter and 
temperature regimesl ^^ * ^^ * ^^ * ^^ For example, high-frequency noise leads to 
the notion of an effective temperature T* of the detecting junction, which 
often exceeds the superconducting transition temperature 

In this work, we first characterize in general the infiuence of current 
fluctuations on a Josephson junction initially in the regime of macroscopic 
quantum tunneling. Consequently, we discuss the sensitivity of a JJ when 
detecting the higher moments of high-frequency noise using switching mea- 
surements. This leads us to investigate the infiuence of weaker low- frequency 
fluctuations that can be considered adiabatic in view of the quantum dynam- 
ics of the detector. Since achieving suitable sensitivity to fluctuations in this 
certain frequency range depends strongly on the electrical environment of the 
detector, we discuss the requirements for a feasible experimental detection 
scheme such as increasing the plasma frequency Up of the detector junction 
or including flltering to limit the accessible noise bandwidth to the desired 
range. 



2. INFLUENCE OF CURRENT FLUCTUATIONS ON A 
JOSEPHSON JUNCTION DETECTOR IN DIFFERENT 
FREQUENCY REGIMES 

In this work we concentrate on discussing how a current-biased hys- 
teretic Josephson junction in the regime of macroscopic quantum tunneling 
can be used to characterize weak fluctuations in its bias current. Such a junc- 
tion can act as a sensitive on-chip detector of current noise since the tunnel- 
ing rate of the phase depends exponentially on current fluctuations. Under 
certain conditions this sensitivity may allow to distinguish non-Gaussian 
features of the fluctuations, which forms a central part of the discussion to 
follow. 

In a general situation when the detector junction is part of an elec- 
trical circuit, these fluctuations of bias current may be of several origins. 
First, equilibrium noise is always present in the circuit, even when external 
current and voltage sources are turned off. For a general linear electrical 
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circuit, the spectral density Si{uj) = dt(/(t)I(0)) exp(ia;t) of these nor- 
mally distributed equilibrium fluctuations is determined by the fluctuation- 
dissipation theorem. As in the previous formula, assuming the distribution 
of the fluctuations to be stationary, the noise power at the detector junction 



where T is the temperature and Y{uj) is the frequency-dependent admittance 
of the electrical circuit seen from the detector. At low frequencies or high 
temperatures with fiw <^ k^T, this reduces to the familiar expression for 
white thermal noise, /S'|°^(w) = ^Gk-^T with G = Re[y(u; 0)] denoting 
an effective conductance. Likewise, at the opposite limit we recover the 
expression for high-frequency quantum noise. 

Besides these equilibrium fluctuations, generally non-Gaussian nonequi- 
librium noise may be present in the circuit as well. For example, in a circuit 
containing a voltage-biased tunnel junction in the normal state, nonequi- 
librium shot noise arises from the stochastic tunneling of discrete charges 
through the junction. 

When both equilibrium and nonequilibrium fluctuations affect the de- 
tector, the total noise power at the detector junction is obtained as the inco- 
herent sum of the different contributions: Si{uj) = Sf^{uj) + Sf°\uj). The 
noise power related to the second cumulant of the shot noise at a finite tem- 
perature in the low-frequency limit is given 

bjff35|hot ^ 2eF2/NC0th(e1//2fcBT), 
where /n = V/Rt<} is the average current through the junction with resistance 
Rf^ at the bias voltage V. Furthermore, F2 is the Fano factor of the second 
moment, relating the microscopic transport properties of the conductor to 
the measurable noise. For a normal tunnel junction with low barrier trans- 
parency, F2 = 1. Similarly, the third cumulant of these nonequilibrium fluc- 
tuations has in the same limit the non-zero theoretical value C3 = F^e^I^^, 
where is the Fano- factor of the third moment .'^^^^'^ The third cumulant 
describes the first-order deviations or asymmetry as compared to a Gaus- 
sian distribution, being therefore one of the first quantities to measure when 
investigating the noise properties of any mesoscopic conductor beyond the 
variance of current. 



The above expressions or their extensions give the spectral densities of 
fiuctuations at an arbitrary frequency and temperature. The response of the 



at frequency uj is given 




(1) 



2.1. Dynamics of a Josephson junction 
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detecting junction to current noise at different frequencies can be explained 
by describing its dynamics using the common resistively and capacitively 
shunted junction (RCS J) -modelP^'^^'^ To be more specific, let us consider 
a JJ with critical current Ic and capacitance C under the influence of an 
external bias current I{t) = /q + 51 {t). Here /q is the time-independent 
average bias current and 51{t) describes the fluctuations. Classically, for 
I < Ic the JJ stays in the supercurrent state and the voltage over the 
junction is zero. On the other hand, a bias current I > 1^ causes it to 
switch to the resistive state, leading to a finite voltage of at least twice the 
superconducting energy gap A to develop over the junction. 

In the RCSJ-model, such a JJ is described in terms of a parallel com- 
bination of a capacitor with capacitance C, an ideal tunnel junction with 
critical current Ic and by a shunt resistance R. The current-phase relation 
of the tunnel junction follows the Josephson relation / = Icsinc^, where 
if is the phase difference of the superconducting order parameter over the 
junction, related to the voltage hy V = hip/2e. 

For most of the discussion to follow, the electrical environment of the 
junction can be described by the effective circuit illustrated in Fig. [H Here 
the current source les contains both a constant bias component /q and a 
time-dependent part related to whose statistical properties are deter- 

mined by the shot noise source. 



3ff 



y(c.) 



C 



Figure 1. An effective circuit model for a J J used as a noise detector. 
Nonequilibrium fluctuations are included in the properties of the effective 
current source producing a current Jeff. Electrical environment of the junc- 
tion is described by the admittance Yiuj) responsible for equilibrium noise. 



For small low-capacitance junctions the charging energy Eq = e^/2C is 
not negligible, leading to a quantum mechanical description of the junctionpSl 
as the charge q and phase over the junction are conjugate variables satis- 
fying q\ = 2ie. In the limit of low dissipation, the junction is described 
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by the RCSJ-Hamiltonian 



cos (f H — — c/p 



(2) 



where E'j = hlc/2e is the Josephson couphng energy. With the above com- 
mutation relation in mind, the behavior of a junction with Eq <^ Ej can 
therefore be visuahzed as that of an imaginary quantum particle with mass 
m = {h/2e)'^C moving in a tilted cosine potential U{(p). At low tempera- 
tures, the (/^-particle is localized in one of the wells of this 'washboard' po- 
tential, performing oscillations at a local minimum at the plasma frequency 
ojp = ^/SEjEcLfQ / h. Here the parameter ipo is related to the bias current / 
changing the tilt of the potential by ipQ = y/l — (//Ic)^- 

Localization of the phase particle corresponds to the junction being in 
the supercurrent state. The state is, however, metastable, and the phase can 
escape from the potential well by macroscopic quantum tunneling through 
the barrier. To estimate the rate of this tunneling process, for currents 
I < Ic one well-barrier section of the potential U{ip) can be approximated 
by a cubic parabolgP^ 

f/M.2Af/(^)7l-ijl), (3) 

4 \ipoJ V ^foj 

where (p is now measured from a minimum of the potential and AU is the 
height of the barrier given by 

Likewise, the tunneling rate r(I) is obtained as 

This follows from the treatment of tunneling out of a metastable cubic poten- 
tial well in the semiclassical limit together with negligibly low dissipation.!^ 
Such a situation arises for small values of the real part of the effective 
shunting admittance Y{uj), corresponding to a high quality factor Q = 
ujpC/Re[Y{Lo)]. 

Assuming that a constant current I is applied for a time At, the finite 
decay rate r(I) leads to the probability 



P{I) = 1 - exp(-r(/)At) 



(6) 
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for the phase to escape from the potential well. After tunneling, the particle 
starts running down the potential hill. For a hysteretic junction with Q ^ 1, 
this corresponds to the junction actually switching to the finite voltage state 
since the particle will not be trapped again into a local minimum until the 
bias current is lowered below the retrapping current close to zero. Hence, 
the quantum tunneling can be experimentally observed by applying repeated 
current pulses of constant height Iq and length At and recording the number 
of times the junction has switched to the normal state. Different values of Iq 
then yield an escape probability histogram P{Io) to be compared with Eq. 

©• 

At higher temperatures, the MQT process is no longer the dominant 
way of escaping from the metastable well. Instead, the energy levels of the 
well have approximately thermal populations and the particle can escape 
over the potential barrier by thermal activation (TA) characterized by the 
rat^ 

Here AC/ is again the height of the barrier introduced in Eq. (j4]). 



2.2. Effect of noise on escape characteristics 

The above view of the dynamics of a JJ allows us to immediately dis- 
tinguish a few different frequency regimes in terms of the response of the 
detecting junction to bias current fluctuations. To see this, let us suppose 
that the switching rate of the detector junction in the presence of noise is 
experimentally determined using the principle presented in Fig. [2l In more 
detail, we assume that constant bias current pulses of height Iq and length 
At injected in a time At/^ are used to obtain each single point on an escape 
probability histogram P{Io). The number N is determined by the desired 
limit on the statistical error in the measurement. Since the total bias current 
I(t) = Iq + 61{t) now contains a fluctuating part, the result is an average 
escape histogram differing from the ideal curve measured without the fluc- 
tuations 51 (t). Then, depending on the frequency content of the fluctuating 
current 61 (t), the pure MQT or TA rate [cf. Eqs. ([5]) and d?])] is modified as 
illustrated by the guidelines of Fig. [3l showing the typical behavior of the 
junction and the relevant characteristic frequencies. 

For noise with frequencies well between AtJ^ and At~^, the detector 
junction is effectively biased by a constant current Iqj during the j:th pulse. 
Consequently, the escape probability is approximately determined as an av- 
erage over the distribution of Iqj probed by the detector when injecting 
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Figure 2. A general measurement scheme to detect current noise witli a J J, 
showing the ideal flow of currents and a micrograph of a typical Al-A10x~ 
Al JJ manufactured using electron beam lithography and 2-angle shadow 
evaporation. In a typical measurement, an average current In is first injected 
through the noise source (a mesoscopic scatterer), whose temporal current 
can then be written as In + SI{t). The DC component of this current is 
returned through the upper arm in the figure. In practise, this is achieved 
by injecting a balancing current —In through the upper line. As a result, 
ideally only the fluctuating current 61{t) flows through the detector junction 
located between points A and B in the figure. This balance of currents 
is verified by monitoring the average current through the detector. For a 
suitably designed circuit, the fluctuations 51{t) pass mainly through the 
detector and do not leak back through the biasing lines. Finally, the escape 
rate of the detector JJ in the presence of the fluctuating current 51{t) is 
measured by injecting N pulses of constant height Iq and length At through 
a third line (middle arm in the figure). The /o~pulses have long leading and 
trailing edges to ensure that the detector responds adiabatically to them. As 
a result of this three-fold current injection scheme, the detector is effectively 
biased by a current Iq + dl{t). The switching of the detector junction out 
of the supercurrent state is detected as voltage pulses between points A and 
B. The effects of non-zero higher odd moments of the distribution of dl{t) 
can be probed by inverting the currents —In and In -|- SI{t), as indicated 
by the alternative signs in the figure. Inverting the average current through 
the noise source changes the sign of the third cumulant of the fiuctuations, 
resulting in differing average escape rates of the detector. 

N current pulses to measure a single point on a histogram. The effects of 
noise at these low frequencies are considered for example in Ref. [25l and 
will not be discussed further in this work. Above this frequency range, weak 
current fluctuations with frequencies up to the order of cjp belong to the so- 
called adiabatic regime. These quasistationary fluctuations have low enough 
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Figure 3. Typical behavior of a J J in the presence of current fluctuations 
belonging to different frequency ranges as discussed in the text. The fluc- 
tuations 51{t) cause the measured escape histograms to shift toward lower 
values of Iq. In this work we focus on the wide adiabatic regime spanning 
frequencies / satisfying At~^ ^ / Here At~^ is the inverse length of 

a bias current pulse and Wp is the plasma frequency of the detector junction. 
Current fluctuations in this regime slowly vary the tilt of the potential, and 
the escape rate is obtained as an average of the MQT rate over the distri- 
bution of the current values. On the other hand, high frequency noise with 
/ ~ ^ excites the J J from its ground state, leading to thermal activation 
at an effective temperature. 

frequencies so that the noise 61{t) can be considered merely as a parame- 
ter slowly varying the tilt of the washboard potential around its average.'^S 
In this limit the junction remains in its instantaneous metastable quantum 
mechanical ground state that decays via MQT. Since the tunneling rate is 
strongly dependent on current, the junction can act as a detector of fluctu- 
ations. 

Contrary to this weak low frequency noise, stronger fluctuations with 
frequencies comparable to and thus to the level separation in the potential 
well excite the phase particle to higher energy levels, leading to escape over 
the barrier top. Again, the junction is sensitive to the fluctuations but 
this time the noise affects the rate of effective thermal activation instead of 
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quantum tunneling. 

Although detecting the influence of the third cumulant of adiabatic cur- 
rent fluctuations is the primary goal of the measurement scheme considered 
in this work, we will first focus on the opposite, high-frequency limit. This 
is necessary as the power of noise at these high frequencies sets a condition 
for the validity of the discussion in the adiabatic regime. 

3. EFFECTS OF CURRENT FLUCTUATIONS AT 
NON-ADIABATIC FREQUENCIES 

To understand the influence of high-frequency noise on a J J detector, 
let us briefly review the model introduced in Ref. [131 The effects of the 
second moment at high frequencies can be modeled as resonant excitation 
between different energy levels in the nearly harmonic potential well of the 
detector junction. At low temperatures, this process is mainly driven by 
the nonequilibrium shot noise. The transitions take place approximately 
at the plasma frequency of the junction, and the resulting level dynamics 
can be described in terms of an effective temperature T* leading to thermal 
activation over the potential barrier. In the following, we choose to denote 
the excitation and relaxation rates between the j — l:th and j:th level in the 
well separated by energy hujjj-i by Jjj-i and 7j-i,j, respectively. 

To derive expressions for the transition rates, we start from the Hamil- 
tonian of a JJ biased with a fluctuating current I{t) = lQ + 6I{t) and write 
Tl in the general form as Tl(t) = Ti.o+V(t), where TCq is the time-independent 
system Hamiltonian and V{t) is a perturbation of the form gAf{t). Here 
3 is a coupling constant and A is an operator of the system (the detector 
junction). The matrix elements {j\A{t = 0)|j — 1) of ^4 between states |j) 
and \j — 1) at t = are expressed as A^'^~^. Further, f{t) is a fluctuating 
operator commuting with A with a power spectrum Sf{uj). 

We proceed by expanding the unitary time-evolution of a state {ipiit)) 
in the interaction picture up to third order in the perturbation V{t): 

^[^-^ l^dnVin) -^1^1^ dndrsT [V{n)V{T2)] 

1^ 1^ l^dndT2dTsT[V{n)V{T2)V{Ts)]^ IV'/(0)), (8) 

where the time-ordering operator T is needed as the operators V{t) at differ- 
ent times do not necessarily commute. To study the influence of the second 
moment of fluctuations, we neglect temporarily the last two terms on the 
right hand side of Eq. ([8]). Now, considering the time-evolution of the ini- 
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tial state \tpj{0)) = \j), we first find the probability amplitude aj-ij{t) = 
~ for a transition to state \j — 1). The corresponding ensemble- 

averaged relaxation rate follows then from the standard expressiorP^ 

7.-1. - = fJ|A^'^'-T5;K,_i). (9) 

A similar calculation gives the excitation rate from the j — l:th to the j:th 
level as 

7,.-i = = (10) 

implying that relaxation is mainly determined by the noise power at positive 
and excitation at negative frequencies. 

The detector junctions usually have Ej ^ Eq, and we can write the 
Hamiltonian TCq in the tight-binding approximation of Eq. ([2]) with I = Iq 
and further approximate the potential by a harmonic one. Referring to Eq. 
([2|), we have g = —Ej/Ic, A = ip and f{t) = SI{t). The required matrix 
elements for the harmonic potential are given by j/2fiujpC, yielding Jj-i.j — 
{j /2hujpC)Sj{L0jj-i) for the relaxation and 'Jjj-i — {j /2fiujpC)Si{—ujjj^i) 
for the excitation rate. Here, the total noise power at the detector 

eV 

Si{±iOjj-i) = 2eF2/Aracoth 



2kBT 



2/ia;,,,_iRe[y(u;,, (^coth ± (11) 



is a sum of the contributions from the shot noise source and the equilib- 
rium fluctuations arising from the dissipative circuit surrounding the junc- 
tion. The former is scaled by a factor a: It relates the noise power at the 
source to the noise power at the detector at frequency Wjj-i- The frequency- 
dependent factor a can be calculated by solving the set of Langevin equations 
written for each branch of the circuit, ^"^ but for the present case we may com- 
bine it with F2 into an effective Fano factor F of the complete circuit. With 
this in mind, the transition rates read 

_ jFeIn coth{eV/2kBT) f H.i-i , 

jFeI^coth{eV/2kBT) ( . 

hwpC +^ Q \ 2kBT ^^^^ 



Here, Q = Cujp/Re[Y{ujjj-i)] denotes the quality factor of the detector 
junction. The level dynamics following from Eq. (jl2p is equivalent to that 
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arising from pure equilibrium fluctuations at an effective temperature T* if 
we require 

and ^,.,_,^^^^fcoth^^-lV (13) 
Comparing this expression with Eq. (jl2p . we have, assuming Wjj-i — i^p, 
fceT* ~ ^ ^ ^. (14) 

2arCOth (coth ^ + QFeI.co^(eV/2k^T) ^ 



In the absence of shot noise when In = 0, T* reduces to the actual tem- 
perature T. For bias voltages e\V\ ^ k^T together with hujp » ksT, T* is 
obtained from 

k^T' ~ (15) 

2arc„th(l + ^) 

Consequently, for high noise curr ents with |/n| > hiolC/QFe, Eq. (fTSi) 
gives T* ~ QFe\Ij^\/2kBUJpC = Fe|/N|/2Re[y(cjp)]. On the other hand, if 
the noise currents are still in the high limit but hujp < 2k-QT, the effective 
temperature contains also a term proportional to temperature T: 

k^T^^k^T+^E^. (16) 

A similarly defined effective temperature is employed also in Ref. [161 

Now, if T* exceeds the crossover temperature for quantum tunneling, 
i.e., T* > To = fiujp/27rkB, the decay of the metastable state occurs primarily 
via thermal activation over the potential barrier.l^ Assuming a bias current 
Iq close to Ic, the potential can again be approximated by a cubic parabola 
and the escape rate F is obtained from d?]) with T = T* . The switching 
probability of the junction for a bias current pulse of height /q and length 
At is correspondingly given by 

P(/o) = 1 - exp (-F(/o)At) , (17) 

allowing comparison with experimental switching histograms at different val- 
ues of 7n. For example, the Iq values corresponding to 50% escape probabil- 
ity, P{Io) = 0.5, or the difference in Iq between the 10% and 90% points allow 
for a straightforward comparison of the theoretical model with experimental 
data. 
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Besides describing the influence of non-adiabatic fluctuations, the effec- 
tive temperature sets a hmit for the apphcabihty of the adiabatic model. To 
find the noise current |/Ar,o| corresponding to the crossover temperature To 
we require that 

QFe\Ijsi,o\ 

1 H fy^2(J = ^' ^^^1 

by use of which we can estimate that the currents |/n| have to stay consid- 
erably below 

Mn.oI = = , r? = coth TT - 1 ~ 0.0037 (19) 

for MQT to be the main escape mechanism. This can be used as a first 
approximation of the limit of validity for an adiabatic description of the bias 
current fluctuations. Imposing such a limit is necessary since the effects of 
the third and higher moments are masked by the second-moment induced 
high frequency effects. To keep the currents In in an experimentally rea- 
sonable regime, the factor QF = F2Q(a'p)Cu;p/Re[y(a;p)] in Eq. ([T9|) can 
be adjusted by using a desired flltering circuit, as will be discussed in Sec. 
\6] As an example, let us consider a detector with Ic = 5 fiA and C = 40 
fF which could well correspond to the parameters of the junction in Fig. [2l 
If the junction is biased at Iq = 0.5/c and the limiting noise current /n,o is 
required to be 5 fiA for typical values Q — 10 and F2 = 1, the factor a{ujp) 
related to the flltering should be of the order of 10~^. 



3.1. Higher order eflfects in the non-adiabatic frequency regime 

When deriving expressions for the above relaxation and excitation rates, 
we assumed a harmonic potential. The corrections to the transition rates 
arising from the anharmonicity of the cubic potential can be handled in 
perturbation theory. Moreover, the influence of higher cumulants of current 
fluctuations is revealed by taking into account more terms in the expansion of 
Eq. ([8]). Yet, in the previous switching measurements^^ these non-Gaussian 
features were masked by the effects of the second moment, well described 
by the above model with an effective temperature. Other experimental de- 
tection schemes considering JJs as detectors of high frequency fluctuations 
may thus become more suitable in the non-adiabatic frequency regime when 
focusing on the higher order effects. 

The weakness of the third-order effects is evident if we consider all the 
terms in Eq. ([8|) and in addition account for the anharmonicity of the cubic 
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potential in first-order perturbation theory. Using the scattering matrix- 
based in-out orderinji^ of the current operators, the averaged time-ordered 
products of the operators can be expressed in terms of three-current correla- 
tors or their spectral functions analogous to the noise power. Assuming these 
functions to be frequency-independent for the range of relevant frequencies, 
we find in the zero-temperature limit that the shot-noise induced transition 
rate from the initial j:th state to the final j + l:th state is given by 

.... - U - ^)^. PO) 



This is nothing but the first term in Eq. (jl2p . arising from a second order 
contribution in a harmonic potential — the third order contribution vanishes 
identically. On the other hand, an initial superposition state like Kj\j) + 
+ 1) leads to oscillating populations of the different energy levels: 
The frequency of these oscillations corresponds to the level separation woi, 
whereas their amplitude is proportional to the third moment of the noise 
source. For a noise source with Fano factor F^, a part of this third order 
contribution originates from a term (j + 1\T [V {ti)V {t2)V {ts)] This 
contributes to the occupation probability of the j + l:th level by 

r f f dri dT2 drg (j + l|r [V{Ti)V{T2)V{n)] I j) 

Jo Jo Jo 

(j _L 1)3/2 _ 

= ^se^^l ' j\^J^^ In [1 - COS (c^oit)] , (21) 

where oq = Correspondingly, the average transition rate arising 

from this term vanishes, predicting that there should be no detectable net 
signal from the third order resonant transitions in the switching measurement 
scheme discussed above. 

Considering the other recently proposed approaches to detecting higher- 
order cumulants, AnkerholcP^ has used an effective Fokker-Planck equation 
to describe the influence of weak short-correlated non-Gaussian fluctuations 
on a Josephson junction in the regime of thermal activation. The detection 
is again based on the sensitivity of the switching rate out of the zero- voltage 
state to variations in current, and the principle of Fig. [2] can be used to 
reveal the presence of higher odd cumulants. 

Another approach has been to consider a noise source capacitively cou- 
pled to a small JJ in the regime of incoherent Cooper pair tunneling;^ The 
third cumulant of high frequency current fluctuations should then be de- 
tectable by comparing the I-V -curves of the detector junction at different 
values of the noise current. On the other hand, in Ref. [28] the authors 
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investigate the influence of third moment of current fluctuations on a two 
level system. This detector can be experimentally realized as a hysteretic JJ 
SQUID. The master equation approach considers the time-evolution of the 
reduced density matrix of the two level system and predicts third-moment 
induced coherent oscillations between different states, observable for exam- 
ple by studying Rabi oscillations of the system. This result is consistent 
with Eq. ()2ip . bearing in mind that Eq. (|2ip is obtained for a multilevel 
system, not a qubit. A similar master equation approach was employed also 
in Ref. [29|to calculate transition rates induced by the third cumulant. 

Furthermore, a measurement with a mesoscopic conductor parallel to a 
current-biased JJ has been predicted to reveal the fourth cumulant of current 
fluctuations. 1^ This is based on a modiflcation of the rate of macroscopic 
quantum tunneling by the mesoscopic conductor. The principle is somewhat 
similar to detecting the influence of the third cumulant in the regime of 
adiabatic fluctuations, as we will discuss below in more detail. 



4. INFLUENCE OF THE THIRD CUMULANT IN THE 
ADIABATIC FREQUENCY REGIME 

The above model with the effective temperature relied heavily on exci- 
tation and relaxation between different energy levels in the potential well of 
the detector junction. In this section we concentrate on fluctuations of lower 
frequencies that do not excite the detector to higher levels. 

To analyze the adiabatic fluctuations quantitatively, let us flrst recall 
that the current distribution p{SI) around the average current Iq is obtained 
as the Fourier transform p{51) = dke~^^^^ ct)si{k) of the characteristic 

function (j)si{k). With c„ denoting the cumulants of the current, the char- 
acteristic function can be expressed as (j)si{k) = exp(^J^2 Since 
the fluctuations are centered around the average current, (61) = ci = 0. 
Consequently, the second and third cumulants are given by C2 = {61 (t)'^) 
and C3 = {61{t)^), and for simplicity we assume the distribution to be sta- 
tionary. Further, assuming the fourth and higher cumulants to be negligibly 
small, we concentrate on the influence of the third cumulant. Truncating 
the sum in the characteristic function at n = 3, one flnds the corresponding 
approximation to the probability density as 

p^SI) ^ -L= (1-^51 + ^5 A exp (-^-^\ . (22) 

This result is valid when the skewness 7 of the current distribution is small, 
i.e., when 7 = c^/c^"^ <C 1.'^ 
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Now, for a J J biased by a current / close to the critical current Ic, the 
MQT rate r(I) is given by Eq. ([5]). Writing the current as I{t) = Iq + SI{t), 
the escape probability for a current pulse of height Iq and length At is 
obtained from 



1 — exp 



At 



r(/o + si{t))dt 



(23) 



For adiabatic current fluctuations 61{t) with frequencies above the inverse 
pulse length At^^ but well below the plasma frequency u;p/27r, we can assume 
the fluctuations to be ergodic. In this case the time average of the escape rate 
can be replaced by an ensemble average over the distribution of fluctuations, 
giving (r) as 

/•OD 

(r) ~ / d5Ip{5I)r{Io + 5I). (24) 



(25) 



This corresponds further to an escape probability 

P(/o) = l-exp(-(r)At), 



which can be directly compared with experimental escape histograms. 

For data fitting we obtain values of (F) by numerical integration but for 
illustration we can write a second order approximation in 51 for F(/o + ^-^) = 
exp [ln^(/o + SI) - B{Io + 61)], giving together with Eqs. ([22]) and ^ the 
result 



(r) 



exp 




C3_A 

2clB 



+ 



_C3_ 

64 



+ 



C3 



A 



24 52 



(26) 



Here Fq is the unperturbed tunneling rate given by Eq. ([5]) with I = Iq and 
we introduced the parameters 



dl 



lnAiI)-B{I) 



i=io 



and 



B 



1 

C2 



dP 



In A{I) - B{I) 



(27) 



(28) 



I=Io 



This approximation is naturally valid only for small variations of current. 
Moreover, neglecting the current dependence of the prefactor A and taking 
only first order corrections in the exponent B{I) into account, Eq. (|26p 
reduces to 
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Figure 4. An example of the shift of the escape histograms compared to the 
ideal noiseless case when adiabatic fluctuations are present at the detector 
junction. The left group of three curves corresponds to a Gaussian noise 
source with C3 = and a non-Gaussian source with equal C2 but C3 7^ 0. 
Additionally, as discussed below in more detail, the arrow indicates the shift 
A/ in bias currents corresponding to a fixed threshold probability (50%) for 
opposite signs of C3. The detector is assumed to have /c = 1 /iA and C = 20 
fF, whereas C2 and C3 are calculated for a current /n = 0.5 /iA in an ideal 
circuit with a cutoff frequency at 0.5 Wp. 




When the approximation of Eq. ()22p is valid, we see that the main effect of 
a nonzero current J/v through the shot noise source is to change the average 
tunneling rate, resulting in a shift of the escape probability histogram. This 
is illustrated in Fig. H] for a Gaussian source and two non-Gaussian sources 
with third cumulants of opposite signs. The current corresponding to a fixed 
switching probability is clearly reduced compared to the noiseless histogram. 
In Fig. [5] such current suppression is shown for the current corresponding to 
50% switching probability when In is varied. 

For different polarities of /at, the change in (F) has the same magnitude 
but different direction. In first approximation, this effect is caused solely by 
the third cumulant. For numerical evaluation of the escape rate asymmetry 
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Figure 5. The bias current Iq corresponding to a 50% switching probabihty 
as a function of In, the average current through the noise source with a 
third cumulant C3 = within the adiabatic model. The detector is assumed 
to have Ic = 2 /xA and C = 20 fF. Furthermore, the circuit around the JJ 
is assumed to have a flat frequency response up to a cutoff at ra;p/27r, as 
shown in the inset. The curves from top to bottom correspond to different 
values of r between 0.5 and 0.8 in steps of 0.1. 




between different polarities of In-, we obtain 

(r)± = ^= / dxTo{h + V^x) exp(-xV2) 
V 27r J -00 

/oo 
dxFoClo + ^/c^x)x(l - xV3) exp(-xV2), (30) 
-00 



±1 ' 



which follows from Eqs. (|22p and (j24p . This can then be used to evaluate 
the relative asymmetry of (F) between different polarities of /at, defined by 
the expression 

AF _ (F)+ - (F)_ 



F,.e i((F)+ + (F)_)' 



(31) 



The above result can further be used when approximating how a certain 
point on the histogram corresponding to a fixed switching probability moves 
as the sign of /n is changed. This shift A/ of a threshold current Iq is 
illustrated in Fig. [6l The value of A/ is numerically easily obtained, e.g., 
by using a bisection method. In the linear approximation, this shift due to 
the third moment is determined by the combined effect of the asymmetry 
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Figure 6. The shift AI/Ic of a threshold current as a function of |/n|, 
the magnitude of the average current through the noise source. AI is the 
difference of the bias currents Iq corresponding to an escape probabihty 
P = X when the polarity of In is changed. The different curves from top to 
bottom correspond to values of x between 0.5 and 0.8 in steps of 0.1. The 
detector is again assumed to have /c = 1 AtA and C = 20 fF. In addition, 
the circuit around the JJ has a flat response up to a cutoff at ruj^j^ix with 
T = 0.5. 

described by Eq. (f3T]) and the slope of the histogram at the average current 
Iq, dP / dI\i=iQ. Thus, for A/ ^ the width of the histogram, we find 

/dP\'^ AT 

A/.(1-P)ln(l-P)(_) _ (32) 

for the shift in the threshold current Iq as Ijy is changed to —In- 

5. EXPERIMENTAL REALIZATION 

In Sec. 12. 2. l and in Fig. [2] we briefly described the principle of detecting 
current fluctuations using a Josephson junction. In this section we aim to 
give a more detailed view, focusing on the electrical circuit model illustrated 
in Fig. El 

The circuit is generally fabricated on a silicon substrate using stan- 
dard electron beam lithography and shadow evaporation of aluminium. For 
measurements the circuit is cooled in a dilution refrigerator down to low 
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Figure 7. Circuit model for the electrical environment of the JJ used as a 
noise detector. The impedances Zd, Zp and are mainly inductive whereas 
Zp and Zb are capacitive. The bias currents Iq and —In are realized by 
applying a bias voltage over a large resistor. When the detecting junction 
switches to the normal state, voltage pulses are measured, e.g., between 
points A and B. 



temperatures (T < 50mK), where MQT is the dominating escape mecha- 
nism in the detector, as desired. The shot noise source is typically another 
small superconducting tunnel junction in the normal state biased by a con- 
stant voltage V. In practise, the noise source is typically realized as another 
Josephson junction with a critical current much smaller than that of the 
detector junction. Referring to Fig. O the average current through this 
mesoscopic scatterer is given by In = V/Rn, where is the resistance of 
the junction at the bias voltage V. For convenience, we separate the DC and 
AC parts of the current /n as In = /n + ^In- Without extra considerations, 
a part of or all of the DC current In will flow through the detector. There- 
fore, in a typical measurement this DC component of the current flowing 
through the scatterer junction is balanced by applying a constant current 
— /n through an inductive line. This current bias is realized by applying 
a voltage over a large resistance at room temperature (i? > 1 MJ7). As a 
result, no DC current flows through the detector junction until a current 
pulse Iq is injected through a similarly prepared line as for the balancing 
current —In- The above procedure allows a more accurate detection of the 
fluctuations for different directions of In, and the balance of DC currents is 
adjusted and monitored for each value of V (or In). 

After the circuit has been balanced, a current pulse Iq is applied. The 
current Iq has the form of a trapezoidal pulse with height Iq. The pulse has 
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long leading and trailing edges to ensure that the detector responds adiabat- 
ically to the changing bias current. Typical pulse lengths vary between 100 
IIS and 10 ms, whereby applying a pulse corresponds to biasing the detector 
with a constant current Iq. Most importantly, the fluctuations caused by the 
shot noise source and possibly attenuated by a low-pass filter pass primarily 
through the detector because of the inductive filtering {Zp and Zb) of the 
lines with the current sources. These fluctuations imposed on top of the 
constant bias Iq can then be probed by sending a large number (typically 
> 1000) of bias pulses of the above kind at a fixed value of Iq repeatedly 
through the detector to produce escape histograms. 

The filter circuit with admittance Yf{uj) in Fig. [7] is essential for guar- 
anteeing the adiabaticity of the current fluctuations, as we will discuss in 
the next section. However, it does not qualitatively affect the measurement 
process described above. On the other hand, the effective admittance Y{u!) 
parallel to the detector junction gives an additional contribution to the tun- 
neling exponent and prefactor of Eq. ([5|). To flrst order, the change of the 
tunneling exponent B can be taken into account for arbitrary Y{uj) following 
Refs. [31] and |32l This results in a correction term AB, which is, however, 
independent of the bias current through the junction. Moreover, the elec- 
trical circuit around the detector affects the cumulants C2 and C3 appearing 
in Eq. ()22p . These are related to the bias current fluctuations over the de- 
tector AId and their distribution, which generally differs from that of the 
fluctuations 5/n at the noise source. 



5L 
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R 



c 



Yj{u) 



Figure 8. The circuit of Fig. [7|at AC frequencies. The impedance Z repre- 
sents the biasing and balancing circuit with the DC sources removed. The 
fluctuating current (5/n can be related to A/d by solving the set of Langevin 
equations as explained in the text. 
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At AC currents, the circuit appears as in Fig. [51 Modeling the noise 
source as a current generator with current 5/n parallel with a resistance 
R^, we have to relate the varying current AId through the detector back 
to J/n- After this, the minimal-correlation estimates for the cumulants can 
be calculated by multiplication and subsequent integration over the relevant 
frequency interval. To achieve this, the current AId is expressed in terms of 

by solving the set of Langevin equations of the circuit!^ For each branch 
j, we first write the total current fluctuation as AIj{uj) = AVj{uj)Yj{uj) + 
6Ij{uj). Here Yj{uj) is the admittance of the single circuit element in branch 
j. This separates the contributions of a fluctuating potential drop AVj{io) 
and those of an actual noise current source, 5Ij{uj). For the circuit of Fig. 
El the only non-zero 6Ij{uj)-teT:m corresponds naturally to (5/n. The set of 
these Langevin equations can then be written as a matrix equation, whose 
solution gives for example the desired relation between A/d and J/n- This 
procedure corresponds to calculating the factor a appearing in the effective 
temperature model discussed in Sec. ET] Furthermore, the relations thus 
obtained can be used to study additional corrections to the third moment 
arising from the electrical circuit and the second moment of fluctuations.'^^ 



6. FILTERING REQUIREMENTS FOR THE 
MEASUREMENT OF ADIABATIC FLUCTUATIONS 

Above we have analyzed the effects on the probability histograms aris- 
ing from the third cumulant of current fluctuations. Adiabaticity of the 
fluctuations is an essential requirement for the theoretical model of Sec. H?] 
to be applicable at all. Low temperature ensures thermal excitation to be 
negligible so that the detector junction stays in its metastable ground state 
in the absence of nonequilibrium fluctuations. On the other hand, in the 
presence of a shot noise source we have to ensure that the electrical circuit 
surrounding the detector junction behaves in a desirable way at high fre- 
quencies. If the noise power is considerable at frequencies near ujp we cannot 
apply the adiabatic model. 

To use the adiabatic model over a wider range of /n, a filtering cir- 
cuit has to be engineered to suppress the high frequency nonequilibrium 
fluctuations. This is illustrated in Fig. [71 by the filter described with a 
frequency-dependent admittance Yf{uj). Ideally, the frequency response of 
the filter and the rest of the circuit surrounding the detector, i.e., the total 
effective admittance Y{uj), should resemble that of a low-pass filter with a 
sharp cut-off below ujp. Due to the admittance Y^lo), the spectral density 
of fluctuations at the detector should behave ideally as in Fig. [9l Since LOp 
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lies generally in the range of 10-100 GHz and the filter circuit has to cover 
effectively a large bandwidth, the high frequencies require us to model the 
circuit using techniques of microwave engineering instead of relying merely 
on a lumped element analysis. 



- 



Figure 9. Idealized frequency responses of the on-chip filtering circuit in case 
of a low-pass filter (solid) or a band-stop filter (dashed). The vertical scale 
corresponds to the squared ratio of the current fluctuation AId^lo) at the 
detector junction to the current fluctuation 5I'^{lo) at the noise source, as 
illustrated in Fig. [8l The characteristics are adjusted so that fluctuations at 
the plasma frequency are suppressed by a desired amount. 



In practise, a band reject fllter with a notch around Wp is easier to 
implement than a good low-pass filter. This frequency response will suppress 
partially the transitions caused by the second moment provided the rejection 
band Auj is wider than the plasma-resonance, Aa; > LOp/Q. Moreover, if the 
detector junction is replaced by a low-inductance DC SQUID with a tunable 
plasma frequency, performing noise spectroscopy of the shot noise source or 
characterization of the fllter circuit becomes possible. 

Considering the experimental realization of such a filter, we manufac- 
tured a sample containing a JJ acting as the noise source and a JJ detector 
separated from the source by a band-stop filter with a stopband centered 
around 30-40 GHz. The theoretical frequency response of the circuit is shown 
as the solid line in Fig. [10] b), whereas Fig. [10] a) shows experimental results 
of the current /q.s corresponding to 50% switching probability as a function 
of |/n| measured using the principle of Fig. [2] The detector is a DC SQUID 
with two parallel junctions with a maximum critical current of 3.7 and a 
total capacitance of 100 fF. Likewise, a JJ with critical current /c,scat — 

1.6 

nA acts as the noise source. The different curves correspond to different 
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Figure 10. a) Measured I0.5 values (open circles) as a function of |In| to- 
gether with theoretical curves calculated using the effective temperature 
modeP^ with QF as a fitting parameter. Different curves correspond to 
different values of magnetic flux through the SQUID loop. The experimen- 
tal points below 1 /iA lie almost at a constant level because the noise source 
is in the superconducting state, b) The QF-va\ues obtained from the fits in 
a) as a function of LOp (right scale). The factor QF is partially determined 
by the frequency response of the on-chip filtering circuit, the theoretical be- 
havior of which is shown as the solid line (left scale). The inset shows a 
schematic of the sample (cf. Fig. [2]), in which the junction acting as a noise 
source and the detector SQUID are separated by a band-stop filter. 



values of magnetic flux through the SQUID loop, and therefore to different 
effective values of Iq and cOp of the detector. For this particular sample, 
the quality factor was very low, Q ~ 2, whereby the energy levels are not 
well separated. A relation between /q.s and |/n| can be obtained using Eqs. 
d?]) and dS]) together with T* from ([ISj), as discussed in Ref. 131 With the 
factor QF as a fitting parameter, one obtains the solid lines in Fig. [TO] a). 
Consequently, Fig. [10] b) shows the variation of these QF~values as tOp is 
altered by changing the magnetic flux. The measured QF~values are gen- 
erally of the order of 0.01-0.1, showing a suppression of almost two orders 
of magnitude compared to typical values in a similar circuit without special 
flltering.J^ This measurement shows that a SQUID-detector can be used for 
noise spectroscopy of the circuit and the noise source. Results of the present 
measurements are, however, not yet well characterized. Improved future 
designs and detectors with higher quality factors should allow to apply the 
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7. DISCUSSION AND CONCLUSIONS 

The rather strict Umits imposed on the adiabatic model can be relaxed 
by using a more general formalism. To treat the nonequilibrium fluctuations 
in the setup of Fig. [1] without necessarily requiring adiabaticity, one could 
approach the problem as quantum tunneling in real time.!^ This corresponds 
to describing the electrical environment and the fluctuations using an ensem- 
ble averaged influence functional and determining the time-evolution of the 
reduced density matrix of the detector junction. Effects arising from the 
third and higher moments of fluctuations can then be analyzed by expand- 
ing the influence functional in terms of the cumulants of current. On the 
other hand, the adiabatic or resonant excitation models considered in this 
work do not take into account the back-action of the detector on the noise 
source: Instead, the nonequilibrium current fluctuations are considered as if 
they were caused by an external, independent driving force. This problem 
can be approached by describing the shot noise source and the detector as 
a coupled quantum mechanical system in terms of an effective action, which 
can further be related to the rate of macroscopic quantum tunneling.UHll 

To summarize, we have analyzed underdamped Josephson junctions as 
detectors of current noise. High-frequency fluctuations generated by a shot 
noise source are described by a thermal activation model with an elevated 
effective temperature. If the junction stays in its ground state in the ab- 
sence of transitions induced by high-frequency noise, the rate of macroscopic 
quantum tunneling is sensitive to the higher moments of fluctuations, too. 
We have discussed the average escape rate for adiabatic fluctuations which 
can be directly compared with measurements. The considered experimental 
scheme allows the detection of non-zero higher odd moments using standard 
switching measurements provided the high frequencies are efficiently filtered. 
Relating the measured tunneling rate asymmetries back to the properties of 
the noise source is a particularly important task, as this involves an accu- 
rate characterization of the frequency dependence of the electrical circuit 
attached to the detector junction. 
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